BI-TP 2005/30 



The impact of QCD plasma instabilities 
on bottom-up thermalization 

Dietrich Bodeker 1 
Fakultat fur Physik, Universitdt Bielefeld, 33615 Bielefeld, Germany 



Abstract 

QCD plasma instabilities, caused by an anisotropic momentum distributions of the 
particles in the plasma, are likely to play an important role in thermalization in heavy 
ion collisions. We consider plasmas with two different components of particles, one 
strongly anisotropic and one isotropic or nearly isotropic. The isotropic component 
does not eliminate instabilities but it decreases their growth rates. We investigate the 
impact of plasma instabilities on the first stage of the "bottom-up" thermalization 
scenario in which such a two-component plasma emerges, and find that even in the 
case of non-abelian saturation instabilities qualitatively change the bottom-up picture. 
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1 Introduction 



It is an interesting and still open question whether in heavy ion collisions it is possible 
to create a medium made of quarks and gluons which is locally in thermal equilibrium, 
which would allow to study thermodynamics of QCD. On the theoretical side it would 
be desirable to have an answer to this question in a limit when one has perturbative 
control, which could be the case for very large collision energies. 

This question can probably be divided into two. The first concerns the earliest 
stage of the collision when partons are "freed" from the colliding nuclei. The second 
is how this system subsequently evolves in time, and whether and when it approaches 
local thermal equilibrium. It should be possible to address the second question in 
perturbation theory, provided the produced partons have sufficiently large energies 
and are sufficiently dense so that multiple scatterings occur. 

In the saturation scenario for high energy nucleus-nucleus scattering pP the freed 
partons are mostly gluons with transverse momenta of order of the so called saturation 
scale Q Aqcd- The distribution function in phase space, or occupation number 
f(x,p) of these gluons is of order 1/g 2 , where g is the QCD gauge coupling (a s = 
g 2 /(4:iT)). The scale that sets the value of g is expected to be Q, so that for large 
collision energies the occupation number is large. Since the typical occupation number 
in thermal equilibrium is of order one, this system is far from equilibrium. 

It is customary to assume that / does not depend on the coordinates transverse to the 
beam direction, which we choose as the z-axis, and that / is invariant under rotations 
around z. This should be a good approximation for sufficiently large nuclei and for 
sufficiently central (head-on) collisions. Furthermore, it is assumed that the gluon 
distribution function is invariant under boosts in z-direction. With these assumptions 
/ only depends on r = \Jt 2 — z 2 , on the magnitude of the momentum perpendicular 
to z, and on the difference of the momentum- and coordinate-space rapidities. 

After the gluons are freed around Qr ~ 1, the system expands and the occupation 
number of the originally produced, so called "hard" gluons drops below 1/g 2 . Without 
production of additional gluons their number density decreases like 1/r. It is not at 
all obvious IH] that such a system thermalizes which would require that the relevant 
interaction rate is larger than the expansion rate 1/r. Consider for instance the rate 
r for 2 — > 2 large angle scattering of hard gluons for which the cross section a is of 
order g A /Q 2 . Taking into account Bose enhancement and using ~ n hjinitial /(<5T) ~ 
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Q 3 / i.9 2 Q T ) one can estimate 

r ~ cm h (l + f h ) ~ g 2 (l + /h)/r (1) 

which is always small compared to the expansion rate. 

In Ref. |1] it was argued that inelastic processes play an essential role for thermaliza- 
tion. In the detailed thermalization scenario of Ref. "soft" gluons are produced by 
bremsstrahlung off the hard gluons. Initially most of the soft gluons have the smallest 
possible energy, but then their momentum rapidly increases in multiple elastic scatter- 
ings with the hard gluons. The soft gluons thermalize among themselves, and are then 
heated by the hard gluons which eventually loose all their energy to the soft-gluon heat 
bath ("bottom-up thermalization"). 

While there is no minimal momentum transfer for massless gluon exchange in the 
vacuum, a medium usually cuts off long range interactions. For particles with nearly 
isotropic momentum distributions this can be accounted for by including the polar- 
ization tensor in the propagators of the exchanged gluons 2 . However, due to the 
longitudinal expansion the typical longitudinal momentum of hard gluons p^ is much 
smaller than the transverse momentum which for a long time remains order Q, so that 
the momentum distribution of the hard gluons is strongly anisotropic. Anisotropic 
distributions can cause so called plasma 3 instabilities which means that some long 
wavelength (|fc| <C Q) gauge field modes grow exponentially if their amplitudes are 
sufficiently small. This is a collective phenomenon which is not visible in the kinetic 
equation framework of this problem used in E3 IH El • The effect has been known 
for many years in plasma physics, and it has been argued jB] that QCD plasma insta- 
bilities can speed up equilibration in heavy ion collisions since they tend to make the 
momentum distributions more isotropic. 

The possible relevance of plasma instabilities to the bottom-up scenario was realized 
in Ref. [7j where also the qualitative difference between QED and QCD plasma insta- 
bilities was discussed. At some point the growth of instabilities is stopped by non-linear 
effects. In QED this happens when the amplitude of the unstable modes has become 
so large that they deflect a particle's momenta by a large angle over the distance of one 
wavelength of an unstable mode. This corresponds to gauge field amplitudes A ~ p/e 

2 The situation is more complicated for magneto-static interactions, but at leading order these do 
not play a role in the scattering processes which are relevant to thermalization. 

3 Here "plasma" refers to a system of quarks and gluons which is not necessarily in thermal equi- 
librium, while sometimes the term "quark-gluon-plasma" is used only for thermalized or almost ther- 
malized systems. 
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where p is a typical particle momentum. Then plasma instabilities have a dramatic 
effect, in particular they are able to isotropize the particle's momentum distribution 
within a very short time. In QCD gluons are self-interacting which may change the 
behavior of instabilities completely. The linear approximation already breaks down at 
much smaller amplitudes A ~ k/g where \k\ <C \p\ is a characteristic wave vector of 
A, and a crucial question is whether these non-linearities stop the growth of instabil- 
ities. In Ref. |H] it was suggested that gluon self-interactions may not saturate the 
instabilities because the system can "abelianize" so that the unstable modes can grow 
until they hit the abelian saturation bound A<p/g. In this case the distribution of 
hard gluons would quickly become isotropic. It has been argued in Ref. |J| that this 
is sufficient for a hydrodynamic description to be applicable even if there is no local 
thermal equilibrium. This could solve the puzzle JU] why hydrodynamic calculations 
successfully describe experimental results for elliptic flow provided that they can be 
used from very early times (r ~ 0.6 fm/c) on, while perturbative estimates for ther- 
malization times are substantially larger (r>2.5 fm/c) [9 . In lattice simulations with 
fields only depending on t and z abelianization was indeed observed [TTJ [T^| , but the 
recent 3+1 dimensional simulations of Refs. jTSJE] indicate that instabilities are sat- 
urated by non-abelian interactions which would mean that their effect is less dramatic 
than suggested in Ref. Nevertheless the instabilities lead to larger infrared gluon 
fields than in an isotropic plasma and their role in the evolution of the system remains 
to be understood. 

So far only the effect of the hard, strongly anisotropic gluons on the infrared modes 
has been considered in the literature. In the bottom-up scenario also the produced soft 
gluons, which have isotropic or nearly isotropic momentum distribution, contribute to 
the polarization tensor and can therefore have an influence on plasma instabilities as 
recently emphasized in Ref. j!5j . In Sec.|2]we consider plasma instabilities in a system 
containing two components, one strongly anisotropic and one isotropic. In Sec. El we 
show how in the early stage of the bottom-up scenario plasma instabilities affect the 
momentum distribution of hard gluons even in the case of non-abelian saturation. 

Notation: 4-vectors are denoted by lower case italics, 3-vectors by boldface. The 
metric is "mostly negative", {g^ v ) = diag(l, —1, —1, —1). 
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2 Instabilities in two-component plasmas 



In the presence of instabilities one cannot naively resum the polarization tensor into 
the gluon propagator and then use it to compute scattering matrix elements entering a 
Boltzmann equation. Instead, one has to fully take into account the dynamics of the low 
momentum gauge field modes. Typically they have large occupation numbers and can 
be described by a classical gluon field = A®T a , where T a are the SU(3) generators 
in the fundamental representation. When the occupation numbers of high momentum 
gluons are small compared to 1/g 2 they can be described as weakly interacting particles 
with a definite momentum. The dynamics of this coupled system satisfies the non- 
abelian Vlasov equations [THl HZj 

(VT = 9 J (2) 



D f + l^U^\ = (3) 



where = — ig[A^, ■] is the adjoint representation covariant derivative, F^ u is the 
field strength tensor and t> M = The distribution functions of high momentum 

gluons are encoded in f(x,p) which is a hermitian N x N matrix with N = 3 being 
the number of colors and f a = tr(T a f). For sufficiently small gauge field amplitudes 
one can write / = ft/N + Sf and linearize with respect to A and Sf. In order to 
see whether there are instabilities one can neglect the x-dependence of / provided the 
growth rate and wave vectors of the unstable modes are large compared to the inverse 
of the time and lenght scales on which / varies. The initial value problem can then 
be solved by spatial Fourier transformation and one-sided Fourier transformation with 
respect to time, A(k ) = J °° dte lhot A{t) where the frequency k is in the upper half of 
the complex plane JE] 4 - By eliminating Sf(k) one obtains an equation for A(k) which 
takes the form 

k 2 g» v - k»k v + W v {k^\ A v (k) = ®»(k) (4) 

with the polarization tensor 

2 f d 3 p df(p) ( _ v^t 



W»(k) = 2Ng 2 / 7^-^ v*6 



(27r) 3 dp 1 \ v ■ k 

d 3 p f(p) ( W + W v»v"k* . 



4 This is of course nothing but a Laplace transformation with a different convention for the variable 
conjugate to t. 
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In Eq. (J3J) is an analytic function of k° off the real axis which depends linearly on 



solve Eq. (jlj) for A(k). If f(p) depends on the direction of p, A(k) in general has poles 
on the imaginary fc°-axis which yield exponentially growing solutions for A(t, k). 
Following Ref. [7] we introduce 5 



in order to characterize the size of IP"'. The physical interpretation of m 2 depends 
on further properties of IP". For example, in an isotropic plasma gluon fields with 
~ \k\ <m are dynamically screened. For anisotropic plasmas there are unstable 
gauge field modes which grow exponentially with a growth rate of order m or smaller. 
For a mildly anisotropic plasma only modes with | k\ < m are unstable while for strongly 
anisotropic plasmas there are also unstable modes with ^> m. 

In the bottom-up scenario there are contributions to the polarization tensor from 
different scales for p. For example, in the first stage (1 <C Qt <C g~ 3 ) m 2 receives 
contributions from hard and soft gluons which are parametrically of the same size. 
In this case m 2 is the sum of and m 2 , which are defined as in Eq. @, but with 
the p-integration restricted to hard and soft momenta, respectively. The hard gluons 
are strongly anisotropic, while the soft gluons are close to being isotropic. Therefore 
there might be no instabilities if only the soft gluons were present. For the subsequent 
discussion the magnitude of the hard and soft momenta play no role. These only enter 
the values of and m B , so we can simply think of the high momentum gluons as some 
two-component plasma. 

Naively one might think that if one adds some isotropically distributed particles 
the instabilities might disappear. That this is not necessarily the case can be easily 
understood from the criterion for the occurrence of instabilities of Ref. It states 
that there is an instability associated with a given wave vector k for each negative 
eigenvalue of the matrix lim e ^o[k 2 5^ — k % W + IP J (ze, k)]. For an isotropic distribution 
II y (fc) vanishes in this limit, so adding particles with isotropic distribution does not 
affect the occurrence of plasma instabilities 6 . 

We write f = ft + fs and correspondingly IP 1 ' = H£ v + IP/' and we assume that / s 

5 In [J] m 2 is called m 2 ^, and parametrically it is the same as rn^ of Ref. 0- 

6 It is of course possible that the added particles increase the collision rate of the gluons which 
contribute to / to the extent that the collisionlcss approximation, which was used to obtain this 
criterion, breaks down. Then the instabilities might disappear. 



the initial values at t — for A(k), A(k), and 5f(k,p). After fixing the gauge one can 




(6) 
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only depends on \p\. Then Tig" has exactly the same form as a polarization tensor in 
thermal equilibrium [20] with the Debye mass mo replaced by v2 m s . 

We are interested in hard gluon distributions for which the typical p z is small com- 
pared to its transverse momentum, i.e., the typical ^-component of v in Eqs. (J^J) is 
much smaller than one. Detailed discussions of the resulting polarization tensor can 
be found in Refs. [3 El- There are no instabilities for \k\ ^> /c max with 

fcmax ~ — (7) 

K 

To keep the present discussion as simple as possible we consider <C \ko\ and 
v h\kz\ ^ l^o I so that we can approximate v ■ k ~ k in the second line of Eq. 
We assume that /h is invariant under p — > — p and under rotations around the z-axis. 
Then one finds 

ng'(fc) ~ (<? j - <W Z ) n£ (k) (8) 

with 

n £ (fc) = ^fc§ + *! 



Note that this is negative when k = i^y and 7 2 < fc 2 . With the same approximations 
the soft contribution to IP- 7 can be written as 

rrf (&) ~ (<^' - 5 lz 5 jz ) n, (*) + <w'*nf (jfe) (io) 

Then there are two unstable modes with polarization approximately orthogonal to z 
and their growth rate 7 is determined by 

7 2 + fc 2 + n ± (i 7 ,fc) = (11) 

where II ± = IT^ + IT^ . As we have discussed above, including IT^ in Eq. (jTTJl does not 
eliminate instabilities but it decreases their growth rate because both IT^ (i 7 , k) and 
ri s (ij, k) increase monotonically with increasing 7 and because II S (ry, k) > 0. 

First consider the case m s ~ m-h- For \k\ ~ mh there is only one scale in Eq. ()1 lj) 
and the instabilities grow with a rate of order m^. When k max 3> |fe| 3> m^, one would 
have 7 ~ mh/v2 if the soft gluons were absent. Since the soft gluons decrease the 
growth rate we have 7 <C |fe|. Then one can approximate 

n s X (fc)^-^m s 2 ^ (N«|fe|, lm(fc°)>0) (12) 
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Using Eqs. Q and (fP2j) and neglecting 7 2 <C fc 2 , Eq. (fTTj) becomes a cubic equation 
for 7, 

* > -5"4? + Hli[ = ° (13) 

Since 7 <C |fe| the third term on the left hand side is small compared to the second and 
can be neglected which again gives 

7~— (m h ~ m s , m h < |fe| < A; max ), (14) 
V2 

i.e., in this regime the effect of soft gluons on the growth rate is negligible. 

Now we consider what happens if m s 3> m^. We will find growth rates which are 
small compared to |fe|. Therefore we can use Eq. (jlUjl to compute 7. First consider 
small |fc|. We neglect the first term in Eq. ()13|) which will be justified in a moment to 
obtain 

7-f^y /3 ifci as) 

Now we see that one can indeed neglect the first term in Eq. (|13|) as long as \k\ -C 
(m s /mh) 2//3 r?T.h. In this range the growth rate is reduced by the soft gluons and it 
increases linearly with |fe|. Going to larger \k\ all terms in Eq. (|13|) become of the 
same order of magnitude and the full cubic equation for 7 needs to be solved. Finally, 
when \k\ 3> (m s /mh) 2//3 mh one can neglect the last term on the left hand side of (JT3J) 
which gives 

7 ~ (m h < m s , m h /3 m 2/3 < |fe| < fc max ) (16) 

Thus, in this case the effect of soft gluons on the growth rate is negligible. Note that 
the range in which Eq. (J16|) holds only exists when <C (mh/m s ) 2 / 3 . 



3 Plasma instabilities and bottom-up scenario 

We will now address the question whether and how plasma instabilities affect the 
parametric estimates of the bottom- up scenario [S]. It is clear that they have a dramatic 
impact if they can grow as large as in QED, i.e., if they do not saturate until A ~ p/g 
(H]. In this case the unstable modes grow so large that they change the momenta of 
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the hard particles by an amount of order one within a time of order m" 1 , which would 
lead to a quick isotropization. We will assume that the instabilities already saturate 
when A ~ k/g, as indicated by recent lattice simulations Then the effect of 

the instabilities on the hard gluons is much less dramatic. 

We only discuss the first stage of the bottom-up scenario when the hard and the soft 
contributions to IP^ are of the same order of magnitude, i.e., m s ~ ~ m and the 
distribution of soft gluons is nearly isotropic. We have seen in Sec. 121 that there are 
plasma instabilities irrespective of the numerical factors in m s and m^. 

The typical longitudinal momentum of hard gluons decreases due to the expansion. 
If they were free streaming, one would have p^ oc r _1 . On the other hand, pf; increases 
due to multiple scattering. Combining these two effects one obtains pi oc r" 1 / 3 j^J. We 
will now consider how the instabilities affect this estimate. 

In the linear regime the unstable modes cannot be described as particles, but rather 
as a coherent classical field. When its amplitude becomes of order m/g its equation of 
motion becomes non-linear. Due to the complicated interaction the coherent nature of 
the gluon field will at least partly be lost. Parametrically the occupation number of 
the ~ m modes is then / m ~ 1/g 2 . A quantitative description of their dynamics by 
a Boltzmann equation is not possible because it would require that f m <C l/g 2 - But 
since this is right at the border of validity of the Boltzmann equation one should still 
be able to use it for parametric estimates which is what we will do in the following. 

Consider 2 —>■ 2 scattering of a hard gluon with a ~ m gluon with a momentum 
transfer of order m. The corresponding cross section a is of order g 4 /m 2 . Therefore 
the rate at which a hard gluon experiences such collisions is 

dN co \ 

—, crn m f m ~ m (17) 

GIT 

where we have used n m ~ m 3 /g 2 for the number density of |fc| ~ m gluons and where 
f m entered as a Bose enhancement factor. The hard gluons experience many random 
collisions each with a momentum transfer of order m so that 

iplf - N col m 2 ~ r^m 2 (18) 

In the bottom-up scenario m 2 ~ g 1 n^jQ ~ Q 2 /{Qt) where is the number density 
of hard gluons. Thus these collisions lead to 

pi ~ QiQr)- 1 ^ (19) 
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which is larger than in the bottom-up scenario where ~ Q{Qt)^ 1 ^ . 

We have thus seen that plasma instabilities indeed lead to a certain isotropization 
in an expanding gluon plasma, even if they are saturated by non-abelian interactions. 
But the isotropization is far from complete, and whether the instabilities really lead to 
equilibration is not clear. 

Let us finally use Eq. (JT9*)l to determine fc max of Eq. ((7j), i.e., the largest \k\ for which 
there are plasma instabilities. We have 

< ~ | ~ (Qr)- 1/4 (20) 

which gives 

U ~ (g^(^) V4 ~ Q(Qt)- 1/4 (21) 

It is interesting that the result (J2"T|) coincides with in Eq. |T9|) . 
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